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ABSTRACT 

The Sun and a Cen A and B are the nearest stars to us. Despite the general agreement 
between their models and seismic and non-seismic constraints, there are serious prob- 
lems pertaining to their interior. The good agreement between the sound speed and 
base radius of the convective zone of the Sun and the solar models is broken apart by 
a recent revision in solar chemical composition. For a Cen A and B, however, it is not 
possible to fit models with the same age and chemical composition to all seismic and 
non-seismic observational constraints. At the age deduced from seismic constraints, 
the luminosity ratio (La/ Lb) of the models is significantly lower than the ratio taken 
from the observed luminosities. Enhancement of opacity as a function of temperature 
is one way to restore the agreement between solar models and the Sun, but such an 
enhancement does not alter the situation for a Cen A and B. The reason is that 
models of both components are influenced in a similar manner and consequently the 
luminosity ratio doesn't change much. In the present study, problems pertaining to 
the interior of these three stars with a single expression for opacity enhancement are 
modelled. The opacity enhancement is expressed as a function of density, ionization 
degree of heavy elements (oxygen) , and temperature. According to this expression, for 
improvement of the models the required opacity enhancement for a Cen A and B at 
log(T)= 6.5, for example, is about 7 and 22 per cent, respectively. The enhancement 
takes place in the region in which pressure ionization is effective, and is higher for low- 
mass stars than for high-mass stars. This result seems to be a possible explanation for 
the serious differences between models and observational results of cool stars. 

Key words: The Sun: interior - stars: interior - stars: evolution - stars: individual: 
a Cen 



1 INTRODUCTION 

Interactions of the constituents of our Universe give shape to 
microscopic and macroscopic objects. In that sense, opacity, 
a measure of interaction between electromagnetic waves and 
particles with mass, is very important for stellar interiors. 
The enhancement of opacity to remove discrepancy between 
models of stellar interioris and seismic and non-seismic con- 
straints on stellar structure is an essential task, because 
opacity influences both energy production and transfer pro- 
cesses dramatically. In the present study, such an enhance- 
ment is tested for improvement of solar, a Cen A and B 
models. 

Opacity enhancement as a function of temperature has 
been the subject of many studies. Simon (1982) considered 
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the influence of opacity enhancement on the period ratio of 
Cepheid oscillations. Korzennik & Ulrich (1989) improved 
the agreement between theoretical and observed frequencies 
of oscillation. Recently, Bahcall, Serenelli & Pinsonneault 
(2004) and Christensen-Dalsgaard et al. (2009) discussed 
removal of the disagreement between the sound speeds of 
the Sun (Basu, Chaplin & Christensen-Dalsgaard 1997) and 
current solar models with recently revised solar chemical 
composition (Asplund, Grevesse & Sauval 2005, hereafter, 
AGS2005). Zdravkov & Pamyatnykh considered the effects 
of opacity enhancement on problems pertaining to hybrid 
stars. The objective of the present paper is to develop a 
model removing the disagreements between interior models 
of the Sun, a Cen A and B and all their seismic and non- 
seismic constraints. For these well known stars, the masses 
of which are similar but different, the opacity enhancement 
is considered as a function of several other physical param- 
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eters besides the temperature. Such an enhancement is in 
accordance with the results of recent studies on non-ideal 
effects in the equation of state (EOS: Saumon Chabrier & 
van Horn 1995, hereafter SCVH). 

By using different types of observational techniques, the 
fundamental properties of a Cen A and B are so well known 
that one can find age of the system from these parameters. 
On the other hand asteroseismology, which is the study of 
oscillatory behavior of stars, puts very important constraints 
on the internal structure of stars such as small and large sep- 
arations between oscillation frequencies. The high precision 
of these constraints for a Cen A and B circumscribes the age 
of the system well. However, the ages found using the fun- 
damental properties and the asteroseismic constraints are 
rather different: 8.9 Gy from the former and 5.5-6 Gy from 
the latter ( Eggenberger et al. 2004; Miglio & Montalban 
2005; Yildiz 2007). At the age derived from the asteroseis- 
mic constraints, the ratio of luminosities (La /Lb) for mod- 
els of a Cen A and B (2.87) is significantly smaller than the 
ratio from the observed luminosities (3.05). Simple opacity 
enhancement as a function of temperature, such as applied 
in the previous studies, affects the structure of a Cen A and 
B in a similar manner and cannot solve the problem with 
the ratio of luminosities. However, the physical conditions 
for a given temperature To vary from star to star. For exam- 
ple, the density at To for a low mass star is in general higher 
than that for a high mass star. Therefore, the required opac- 
ity enhancement should also depend on density. 

Opacity in a medium is a function of cross-section and 
therefore of the ionization degree of ions in equilibrium. 
Equilibrium is the case in which ionization and recombi- 
nation processes balance each other. Whereas the rate of 
the ionization process depends upon the number of energetic 
photons and henceforth the temperature of the medium, the 
density, which is important for non-ideal effects such as pres- 
sure ionization, determines the rate of the recombination 
process. Thus, ionization equilibrium is controlled by the 
density and temperature of the medium. As a result, opacity 
enhancement as a function of temperature is not sufficient. 
Hence, we develop a more comprehensive approach for the 
enhancement. 

For an accurate opacity computation, the EOS must 
be computed very precisely. In addition to the ionization 
degree, the excitation level of an ion is important because 
cross-section of an ion also depends on its excitation level. 
In stellar interiors, however, collisional rates for the excita- 
tion process dominate radiative rates (Mihalas and Weibel- 
Mihalas 1999). This is another indication showing the im- 
portance of density in opacity because the collision rate in 
a medium strongly depends on its density. 

For a realistic opacity prediction, atomic data of high 
quality are required. It is known that although the opacity 
tables are based on extensive models and calculation, the 
calculations do not consider all aspects of the opacity for a 
given temperature and density (SCVH) In the present work 
we assume that the existing opacity tables are uncertain 
at certain densities and temperatures that leave room for 
modifications. 

In Yildiz et al. (2006), to remove the deviation between 
model and observed radii (or colors), a mixing- length param- 
eter is considered. However, enhancement in opacity, which 
controls the energy transfer in the radiative part, increases 



the size of a star. In addition to that, the opacity in the 
core is also very important for the energy generation rate 
because the temperature of the core depends on how much 
released energy is permitted to escape. 

The constraints from oscillation frequencies (y n ,i with 
order n and degree V) inferred from asteroseismic observa- 
tions are expressed in terms of small {v n .i — ^n-ij+2) and 
large (Az^(n) = v n> i — v n -i,l) separations between these 
frequencies. The small separation (Christensen-Dalsgaard 
1988) for / = is defined as 

Si>o2(n) = v„fl - v„-i } 2- (1) 

The turning point of oscillations depends mainly on the 
degree of oscillations I. Whereas the oscillations with I — 
sinks down to the center, the turning points of oscillations 
with I — 1 and I — 2 are about r t = 0.057?* and r t = 0.10.R*, 
respectively (Yildiz 2008). As nuclear evolution proceeds, 
the sound speed throughout the core decreases, and the 
oscillations with I — 1 and I = 2 sink deeper. The result 
of this process is a reduction in the small separation with 
time. For the Sun, a Cen A and B, the oscillations with 
degree I = 0, 1, 2, 3 are detected (Basu et al. 2007; Bedding 
et al. 2004; Kjeldsen et al. 2005). The frequencies of these 
oscillations enable us to compute observational values of 
8voi(n) = (u„o - (fcVi-i.i + v n i)/2) (Kjeldsen et al. 2005), 
5vo2 (n) and 5^13 (n). 

The oscillation frequencies are typically computed as- 
suming adiabatic process. However, this approximation is 
valid for the interior but not for near-surface regions. The 
computation of oscillation frequencies in near-surface re- 
gions requires perturbation of the energy equation. Unfortu- 
nately, perturbation of the energy equation results in a large 
uncertainty in the convective flux (Christensen-Dalsgaard et 
al. 1996). Kjeldsen Bedding & Christensen-Dalsgaard (2008) 
derived an empirical correction for the near-surface offset 
between observed and computed oscillation frequencies by 
assuming a power law for the offset. They showed that this 
offset is a function of frequency. Use of small and large sep- 
arations or their ratio (Roxburgh and Vorontsov 2003) min- 
imizes the effect of the near-surface regions. 

The small separations are a function of n. In the present 
study, we use their average values over the observed oscilla- 
tions for each star, for comparison with model values. The 
average values of small (Svoi, S V02/Q and 5^13/ 10) and large 
separations (Auo) for the Sun are calculated using the fre- 
quencies of an order between n=ll and n=27. The solar 
small separations are given in Table 1. The ranges of n 
for a Cen A and B are n= 16-25 and n=20-27, respectively. 
For 1=0, SV02/6 values of a Cen A and B are 1.05±0.12 
and 1.69T0.08, respectively. The small separations of some 
stars are less than zero for some values of n (see e.g. Soriano 
& Vauclair 2008), however, they are always positive for the 
Sun, a Cen A and B. 

The remainder of this paper is organized as follows. In 
Section 2 and 3, we discuss opacity enhancement as a func- 
tion of temperature in order to improve solar models and 
models of a Cen A and B, respectively. Models with en- 
hancement as a function of several physical parameters in- 
cluding temperature are given in Section 4. Section 5 is 
devoted to possible connection between pressure ionization 
and opacity enhancement. Finally, we give concluding re- 
marks in Section 5. 
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Figure 1. Fraction of OV in total number of oxygen ions within 
the Sun (solid line), a Cen A (thin dotted line) and B (thick 
dotted line). 

2 TEMPERATURE-DEPENDENT 

ENHANCEMENT TO IMPROVE SOLAR 
MODELS 

The degree of ionization has a particular importance in the 
resistivity of matter against electromagnetic radiation. In 
the central regions of stars, the temperature is high enough 
for almost complete ionization of hydrogen, helium and the 
most abundant heavy elements. Therefore, in these regions 
the cross-section of ions is very small and the opacity is 
low. High ionization occurs in the central regions of stars 
while low ionization or neutrality takes place in the outer 
regions (see Section 5). Specific ions occupy specific parts 
of the stars. This phenomenon is depicted in Fig. 1. The 
fraction of OV (oxygen ion with four electrons) in the total 
number of oxygen ions is plotted with respect to log(T) for 
the models of Sun (solid line) and a Cen A (thin dotted line) 
and B (thick dotted line). The distribution of OV is likely 
to follow a Gaussian distribution. 

Therefore, we write the artificial enhancement of opac- 
ity as 

tenh = KOPAL + 5n. (2) 

where 5k is always positive and represented by a Gaussian 
function. 

8k = kopalci lo-^m-csf (3) 

In equation (3), kopal is the opacity computed from the 
OPAL tables (Iglesias and Rogers 1996) and ci, C2 and C3 
are free parameters to improve the solar models. c\ is ampli- 
tude of enhancement. While C3 gives temperature (log(T)) 
at which the enhancement is maximum, C2 shows how the 
enhancement is localized. 

Solar models are constructed with metallicity Z — 0.016 
for different values of ci, C2 and C3. The basic properties of 
the models are listed in Table 1. In the first row, the standard 
model (Yildiz 2008) is given (ModST). In the second, third 
and fourth columns, the numerical values of the free param- 



eters ci, C2 and C3 are listed, respectively. For ci = 0.25, 
the solar models for three different values of C2 are given in 
the second (ModKl), third (ModK2) and fourth (ModK3) 
rows of Table 1. A decrease in C2 reduces the base radius 
of the solar convective zone (fbcz), however it augments the 
helium abundance in the solar envelope (Y s ). The fifth row 
gives the ModK4 solar model. For this model, opacity en- 
hancement is more localized to the regions just below the 
convective zone. Its Y s and rb cz are in very good agreement 
with the observed values. In addition to the good agreement 
between the observed and predicted base radius of convec- 
tive zone, the small separations (Svoi, 81/02 and 81/13) and 
the sound speed concur reasonably well. The relative sound- 
speed differences between these models and the Sun (RSSD) 
are plotted in Fig. 2. The RSSD for the standard solar model 
(ModST) is about 1.7 per cent and decreases with enhanced 
opacity below the convective zone. Thin dotted, thick solid, 
and thick dotted lines in Fig. 2 show ModKl, ModK2, and 
ModK3, correspondingly. Among these, ModK2 (02 = 4.5) 
is in better agreement with the Sun than the others. An- 
other solar model (ModK4) in which enhancement of k is 
more localized to below the convective zone (ci = 0.30 and 
C2 = 20) gives better results for Ys and rbcz 

than the solar 

models with ci = 0.25. 

A better model (ModKM in the sixth row of Table 1) 
is constructed with the mean of enhancements of ModK2 
and ModK4. This ModKM model is the best of the solar 
models regarding all seismic constraints and including the 
small separations given in the last three columns of Table 
1. The RSSD for the ModKM model is less than 0.2 per 
cent (solid line with diamonds in Fig. 2). Also shown in Fig. 
2 is the RSSD for the solar model constructed with opacity 
enhancement (dotted line with +) as given in equation (4) of 
Bahcall et al. (2004). The maximum RSSD for the Bahcall 
et al. (2004) model is about 0.6%. Its basic properties are 
given in Table 1 (ModB4). 



3 MODELS OF a CEN A AND B WITH 

ENHANCED OPACITY AS A FUNCTION 
OF TEMPERATURE 

Treating the enhanced opacity as a function of temperature 
improves the solar models. If we apply any of the enhance- 
ments discussed in the previous section, for example that of 
Model 3, to models of a Cen A and B with chemical compo- 
sition given by Feltzing & Gonzalez (2001), the fundamen- 
tal properties (L, R, T e g) of the models of a Cen A and B 
are essentially unchanged. The same is true for the opacity 
enhanced as in Bahcall et al. (2004). The basic properties 
of models of a Cen A and B with the enhanced opacity of 
Bahcall et al. (2004) are listed in the fourth and fifth rows 
of Table 2 (SIS23KBah), respectively. These models are ob- 
tained with the parameters of Model SIS23 of Yildiz (2008). 
For comparison, the basic properties of Model SIS23 are also 
listed in the first and second rows of Table 2. Despite the neg- 
ligibly small differences between the luminosities and radii 
of the models with and without opacity enhancement, the 
small separations of Models SIS23KBah for a Cen A and B 
(with a moderate enhancement of opacity) are slightly dif- 
ferent from those of Models SIS 23. This is consistent with 
the result of previous section that the opacity enhancement 
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Table 1. Basic properties of the solar models with the chemical composition given by AGS2005. The last line is for the observed values. 
The uncertainty in R c is 0.001 Rq (Basu & Antia 1997). The observed values of the surface helium and metal abundances are taken 
from Basu & Antia (1995) and AGS2005, respectively. The unit of small separations Svoi, Suq2, and <5fi3 (Basu et al. 2007) is /iHz. 



Model 


Cl 


C2 


C3 


X 


a 


Pc 


T c 


x c 


Z 3 


Y 3 


Rc/Rq 


Suqi 


<5f02/6 


Svis/W 


ModST 








0.70975 


1.820 


146.78 


15.65 


0.355 


0.0124 


0.244 


0.733 


1.879 


1.686 


1.739 


ModKl 


0.25 


18.0 


6.30 


0.70901 


1.871 


146.85 


15.65 


0.355 


0.0128 


0.248 


0.7162 


1.821 


1.662 


1.716 


ModK2 


0.25 


4.5 


6.30 


0.70506 


1.890 


148.46 


15.70 


0.349 


0.0129 


0.252 


0.7130 


1.811 


1.646 


1.701 


ModK3 


0.25 


3.0 


6.30 


0.70200 


1.906 


149.35 


15.72 


0.346 


0.0129 


0.255 


0.7125 


1.800 


1.638 


1.696 


ModK4 


0.30 


20.0 


6.32 


0.70881 


1.877 


146.93 


15.65 


0.354 


0.0128 


0.248 


0.7130 


1.821 


1.661 


1.715 


ModKM 








0.70684 


1.883 


147.76 


15.68 


0.352 


0.0128 


0.250 


0.7130 


1.814 


1.653 


1.709 


ModB4 








0.70898 


1.857 


146.80 


15.65 


0.355 


0.0128 


0.248 


0.7140 


1.811 


1.688 


1.740 


Obs 


















0.0122 


0.25 


0.713 


1.813 


1.653 


1.683 



Table 2. Model properties of a Cen A and B. <5^o2 and A^o are in units fiHz. The error in Ai/o is 0.1 /iHz. The uncertainty of 5vo2 is 
0.12 ^tHz for a Cen A and 0.08 fMz for a Cen B. 



Star 


L/Lq 


R/Rq 


T c ft 


X 


Zo 


a 


t(W 9 y) 


5^02/6 


Au 


MODEL 


A 
B 


1.549 
0.569 


1.223 
0.863 


5828 
5400 


0.7129 
0.7129 


0.023 
0.023 


1.710 
1.980 


5.90 
5.90 


0.91 
1.70 


106.7 
164.4 


SIS23 
SIS23 


B 


0.507 


0.864 


5246 


0.7129 


0.023 


1.820 


5.90 


1.59 


164.7 


ModK3 


A 
B 


1.535 
0.567 


1.224 
0.863 


5824 
5395 


0.713 
0.713 


0.023 
0.023 


1.71 
1.98 


5.90 
5.90 


0.97 
1.56 


106.4 
164.5 


SIS23KBah 
SIS23KBah 


A 
B 


1.545 
0.506 


1.226 
0.862 


5817 
5247 


0.7066 
0.7066 


0.023 
0.023 


1.881 
1.930 


7.20 
7.20 


0.65 
1.45 


106.4 
165.0 


ModK5z23 
ModK5z23 


A 
B 


1.546 
0.507 


1.226 
0.863 


5817 
5249 


0.7066 
0.7066 


0.023 
0.023 


1.731 
1.791 


5.90 
5.90 


0.97 
1.60 


106.1 
165.0 


ModK6z23 
ModK6z23 


A 
B 


1.545 
0.506 


1.223 
0.863 


5818 
5246 


0.7135 
0.7135 


0.021 
0.021 


1.676 
1.775 


5.56 
5.56 


1.04 
1.64 


106.7 
164.9 


ModK%z2\ 
ModK%z2\ 


A 
B 


1.543 
0.507 


1.224 
0.863 


5819 
5248 


0.7126 
0.7126 


0.019 
0.019 


1.743 
1.874 


5.58 
5.58 


1.01 
1.62 


106.7 
164.9 


ModKSzl9 
ModK8zl9 


A 
B 


1.544 
0.507 


1.224 
0.863 


5824 
5250 










1.05 
1.69 


105.5 
161.5 


obs 
obs 





Table 3. Basic physical conditions at the bottom of the convective zones of the standard solar and a Cen A and B models and amount 
of opacity enhancement there to improve these models. 



Star 


Mc/Mq 


T bcz (MK) 


Pbcz(gcm 3 ) 


Cl 


a Cen B 


0.934 


2.69 


0.571 


0.425 


The Sun 


1.000 


1.97 


0.150 


0.250 


a Cen A 


1.105 


1.73 


0.068 






Table 4. Basic properties of the solar models with enhanced opacity as a function of p, T, and ionization degree of oxygen. See also 
Table 1. 



Model 


X 


a 


Z 


T c 


x c 


Z s 


Y s 


R c /Rq 


Suoi 


<5^02/6 


5v 13 /W 


ModK5 


0.70174 


1.875 


0.0160 


15.71 


0.349 


0.0128 


0.255 


0.7167 


1.801 


1.635 


1.691 


ModK6 


0.70163 


1.885 


0.0153 


15.67 


0.353 


0.0124 


0.257 


0.7111 


1.811 


1.688 


1.740 


ModK7 


0.70106 


1.898 


0.0155 


15.63 


0.357 


0.0125 


0.248 


0.7130 


1.811 


1.660 


1.716 


ModK8 


0.70854 


1.892 


0.0150 


15.61 


0.357 


0.0121 


0.250 


0.7130 


1.809 


1.655 


1.710 


Obs 












0.0122 


0.25 


0.713 


1.813 


1.653 


1.683 
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Figure 2. The relative sound-speed difference between the solar 
models given in the first three rows of Table 1 and the Sun is 
plotted with respect to the relative radius. The thin solid line 
shows the standard solar model with Zq = 0.016 (Yildiz 2008). 
The other models are with enhanced opacity and given in Table 
1 (see the text). 

influences the sound-speed profile, at least in the central 
regions. However, the effects on a Cen A and B are in oppo- 
site directions; the enhancement decreases 5uo2 of a Cen B 
but increases Svo2 of a Cen A. Therefore, the temperature- 
dependent opacity enhancement concept does not improve 
the a Cen A and B models. 



4 THE GENERALIZED ENHANCEMENT OF 
OPACITY 

The agreement between observed and model predicted lu- 
minosity values of a Cen B may be improved by increasing 
the coefficient ci. For this purpose, we set ci =0.425. The 
results from the ci =0.425 model are given in the third row 
of Table 2. Since model SIS23 of a Cen A is already in very 
good agreement with seismic and non-seismic constraints, 
opacity enhancement is not needed. Hence, the opacity en- 
hancement factor is assumed to be zero for a Cen A. Next, 
the ci factor is considered to be a function of stellar mass. 
The basic physical conditions at the bottom of the convec- 
tive zones of the standard solar and a Cen A and B models 
are listed in Table 3. The second column in Table 3 lists for 
the masses of the stars. The third and the fourth columns 
report the base temperature and density of the convective 
zones. In the fifth column, the values of coefficient ci, nearly 
the maximum value of 5k/ k, are given. Despite the very dif- 
ferent chemical composition of the models for the Sun and 
the components of a Cen, the relation between the factor ci 
and stellar mass is nearly linear (see Fig. 3). This relation 
can be attributed to density dependence of ci . The densities 
below the convective zones of the Sun, a Cen A and B are 
very different from each other. For example, the density at 
the base of the convective zone of a Cen B is about four 
and eight times higher than those of the Sun and a Cen A, 



1 1 1 1 1 1 1 1 1 ^ 1 

0.92 0.96 1.0 1.04 1.08 1.12 

M/M 

Figure 3. ci and base density of convective zone as a function 
of stellar mass. 
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Figure 4. The relative sound-speed difference between the solar 
models and the Sun. 



respectively. Since the base temperatures of the convective 
zones of these stars are much closer to each other than the 
densities, ci can be considered as a function of density only. 
The density ranges of the regions in which opacity is effec- 
tively enhanced are very different for the Sun and a Cen A 
and B, and move toward higher densities for low mass stars. 

As stated above, the dependence of ci on stellar mass 
can be related to density in the region below the convective 
zone, which differs from star to star. Therefore, some ex- 
pressions for opacity enhancement as a function of density 
and temperature, which improve the solar and the a Cen A 
and B models, can be tested. For this purpose, ci is plot- 
ted as a function of stellar mass in Fig. 3 (solid line with 
diamonds). Also plotted in Fig. 3 is the base density of the 
convective zones (pbcz) for each star. There is a correlation 
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log(T(K)) 



Figure 5. The enhanced opacity in the Sun with respect to 
log(T). The base temperature of the convective zone (7b cz ) of 
the Sun is given by log Tt, cz =6.34. The maximum enhancement 
{5k/k) for the solar model is about 0.1. 
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Figure 6. The enhanced opacity in a Cen A and B with respect 
to log(T). The base temperatures of the convective zones (Ti, cz ) 
of a Cen A and B are different: log Tt, cz of a Cen A and B arc 6.31 
and 6.50, respectively. The enhancement (Sk/k) at logT| 3cz =6.5 
is about 0.067 for model of a Cen A and 0.22 for model of a Cen 
B. 

between pbcz and ci. We therefore construct a solar model 
with Sk p = 10plO~ 3(logT ~ 5 ' 8) . The characteristics of this 
solar model are given in the first row of Table 4 (ModK5). 
The seismic properties from the new model are not as good 
as those of the solar models with temperature-dependent 
opacity enhancement. However, they are much better than 
those of the standard model ModST. The small RSSD for 
ModK5 is plotted with respect to relative radius in Fig. 4 
(solid line). 



We also construct models for a Cen A and B with 8n p . 
The basic properties of these models (ModK5z23) are listed 
in Table 2. The age from the fundamental properties of a 
Cen A and B is about 7.2 Gyr. At this estimated age, the 
model values of small separations are significantly smaller 
than the observed values. Despite the improvement, the ages 
derived from the fundamental properties and the seismic 
constraints are not in good agreement; the age of the models 
is about 1.5 Gyr greater than the seismic age. 

Sk p given above may solve the inconsistencies pertain- 
ing to the solar interior; however, it does not help with a 
Cen A and B. In order to obtain good agreement between 
the ages of a Cen A and B from fundamental properties 
and asteroseismic properties, the coefficient of Sk p should 
be increased. Unfortunately, such a treatment hurts the 
agreement between the solar model and the Sun. Therefore, 
further parameters are required for opacity enhancement. 

The problem with ci as a linear function of just p is that 
it increases toward the center of stars and becomes maxi- 
mum at the center. On the other hand, at the center of stars 
the opacity is essentially very low as a result of high ioniza- 
tion and therefore no or little enhancement is expected. For 
the radiative interior of cool stars, the effective charge (mean 
ionization degree) of oxygen (z7 e ff,o) is a good indicator for 
opacity. For a completely ionized medium, where Z e ff,o is 
equal to atomic number of oxygen (Zo), a very low value of 
opacity is determined by free-free transitions. In that case, 
we assume no opacity enhancement. Therefore, to minimize 
opacity enhancement near the center, we multiply Sk p by 

(1 f^) 1 ^ 2 ) which is equal to zero for the completely 

ionized medium (see below). 

The EOS routines of the code that we use solve the Saha 
equation (Yildiz and Kiziloglu 1997). The Debye- Hiickel 
approximation is employed for Coulomb effect and EOS of 
Mihalas ct al. (1990, hereafter MHD) is used to compute 
the partition functions. The ionization degree of the ten 
most abundant elements is computed accurately. Oxygen is 
the most abundant heavy element in the Sun and therefore 
opacity is very sensitive to the degree of its ionization, at 
least in the inner radiative region. Therefore, we enhance 
the OPAL opacity as 

^n h = KO PAL(l + 42.5p(l-^) 1 / 2 lO^m-^) 2 ).^) 

Zo 

In terms of the number of oxygen ions with charge i (Ni), 
Z c «,o = X] *-^« ' w here A^o is the total number of oxygen 
ions. 

When it comes to sound speed, the solar model with the 
latest opacity enhancement (ModK6) agrees the best among 
all the solar models. The results from the ModK6 model are 
shown in Fig. 4. The dotted line with diamonds in Fig. 4 
represents this model. The maximum value of the relative 
difference is about 0.15 per cent. The basic properties of the 
solar model ModK6 are given in the second row of Table 4. 
Its Y s , 0V02 and 0V13 are slightly greater than the observed 
values given in the fifth row of Table 4. 

We also construct models for a Cen A and B with 
the enhanced opacity given in equation (4) for Z=0.023 
(ModK6z23) and 0.021 (ModK6z21). The properties of these 
models are reported in Table 2. The small separations (SU02) 
of both models are in agreement with the asteroseismic in- 
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ferences. However, the models with Z=0.021 display better 
agreement than those with Z=0.023. 

We also obtain solar models based on modified version 
of equation (4): 

fl i A c l-l ^eff,Osl/2 1n -17(log(T)-6.32) 2 x ,r\ 
Kenh = KOPAL(l + 4.5p(l ) 10 8 ^ ).(5) 

Zo 

The solar model with this opacity (ModK7) is given in the 
third row of Table 4. Its Y s , 8U02 and 81/13 are in better 
agreement with the observed values than that of ModK6. 

Also given in Table 4 is ModK8, the opacity of which 
is taken as the mean of equations (4) and (5). This model 
is also in very good agreement with the Sun. The models 
of a Cen A and B with this opacity and Z= 0.019 are also 
constructed; ModK8zl9 A and ModK8zl9 B in Table 2 are 
also in good agreement with their seismic and non-seismic 
constraints. 

The opacity enhancement of ModK5s (ModK5 for the 
Sun and ModK5z23 for a Cen A and B) is a function of 
p and T; f{p, T) = Sk p = 10pl0~ 3(log T ~ 5 ' 8) . For ModK6s 
(ModK6 for the Sun, ModK6z23 and ModK6z21 for a Cen 
A and B) and ModK7, the enhancement is as given in equa- 
tions (4)- (5), respectively. The main difference between these 
enhancements is that equations (4)- (5) are also a function of 
ionization degree (effective charge) of heavy elements (oxy- 
gen). Opacity in the (radiative) interiors of the Sun, a Cen 
A and B is essentially determined by bound-free transitions 
(ionization). However, the effect of this process is null if the 
gas is completely ionized. In this case, becomes unity. 

Then, no enhancement is required when the term (1— ) 
in equations (4)- (5) is zero. Since the models with enhance- 
ments with equations (4)- (5) are in good agreement with the 
observed properties of these stars, we deduce that i) simple 
enhancement as a function of p and T is inadequate, and ii) 
much more precise methods for the computation of ioniza- 
tion and cross-sections of heavy elements are required than 
used at present. 

4.1 Comparison of opacities of standard and 

improved models for the Sun, a Cen A and B 

The opacities (log k) of standard and improved solar models 
are plotted with respect to log(T) in Fig. 5. The thin solid 
line shows the standard model (ModST), which has the low- 
est opacity in the outer part of the solar radiative core. The 
thin dotted line and thick solid line represent the improved 
solar models ModK6 and ModK8, respectively. In the core, 
the opacities of solar models are very close to each other; 
in other words, the enhancement is not needed there. In the 
outer parts of the solar radiative core, however, the opac- 
ities of improved models are significantly higher than that 
of the standard model. The opacity difference between the 
standard model and ModK8 is remarkable for log(T) < 6.6. 
The opacity of ModK6, however, is the same as that of the 
standard model only in the central regions (log(T) > 6.9). 

In Fig. 6, we compare opacities of standard (SIS23, thin 
lines) and improved (ModK8zl9, thick lines) models for a 
Cen A and B. The opacities (log/c) in the interior of a Cen 
A (dotted lines) and B (solid lines) are plotted with respect 
to log(T). It should be noticed that the chemical compo- 
sition of standard and improved models are not the same 
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Figure 7. Evolutionary tracks of models for a Cen A and B 
in the HR diagram. While the dotted lines show the standard 
models (SIS23) of a Cen A (thin) and B (thick), the solid lines 
represent the improved models (ModK8) for a Cen A (thin) and 
B (thick). 

and therefore the difference between the opacities is not en- 
tirely due to the enhancement. Due to a lower metallicity 
of improved models than the standard models, despite the 
enhancement the opacities of improved models near the cen- 
tral regions of a Cen A and B are lower than those of their 
standard counterparts. In the outer part of inner radiative 
regions of a Cen A and B, the enhancement is effective for 
log(T) < 6.7 and log(T) = 6.9, respectively. The opacity 
difference between the improved and standard models of a 
Cen B is greater than the difference for a Cen A. This result 
is due to the density and effective charge dependence of the 
enhancement applied to the improved models. 

The bases of the solar and a Cen B models have very 
different densities and temperatures (T bcz ). logT bcz = 6.34 
for the Sun and logTbcz = 6.5 for a Cen B. The enhance- 
ment {Sk/k) at the base of the convective zones needed to 
improve the solar and a Cen B models is about 0.1 and 0.17, 
respectively. 

For assessment of the density dependence of opacity en- 
hancement, it is better to give the amount of enhancement at 
the same temperature. For example, at T = 2.69 x 10 6 K, the 
base temperature of convective zone of a Cen B, the amount 
of opacity enhancement {8k/ k) is 0.067 for a Cen A and 0.22 
for a Cen B (ModK8zl9 with Z = 0.019). For the improve- 
ment of the solar model (ModK8 with Zo — 0.015), the 
amount of enhancement at the same temperature is 0.042. 

4.2 Comparison of evolutionary tracks of standard 
and improved models for a Cen A and B 

Changes in interior opacity significantly affect stellar struc- 
ture as a whole. The opacity enhancement given above al- 
ters evolutionary tracks of cool stars, in particular. In Fig. 
7, the standard and improved models are plotted in a theo- 
retical HR diagram. The observed positions of a Cen A and 
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the effective temperature difference between the enhanced 
and the standard models (T fr, C nh — Toff.st) is -357 K for 0.8 
M and -286 K for 0.9 M . 

We can also compare sizes of standard and enhanced 
models of the same luminosity. For L =0.25Lq, for example, 
the difference between the radii of standard and enhanced 
models of O.8M is about 4.6 per cent. The same amount of 
difference also exists for the 0.9 M model for L=0.41 and 
0.90L Q . The temperature difference is now about 100 K for 
both models for any fixed value of luminosity. As a result, 
opacity enhancement may be a possible explanation for why 
low-mass stars are cooler than (or oversize compared with) 
their interior models (see e.g. Clausen et al. 2009). 
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Figure 8. Evolutionary tracks of 0.8-(thick) and O.9-M0 (thin) 
models with ModK8zl9 values of a Cen B. While the dotted lines 
show the standard models, the solid lines represent the models 
with enhanced opacity. 



B (circle) are also seen. The dotted lines show the standard 
models while the improved models are represented by the 
solid lines. The evolutionary tracks comprise main-sequence 
(MS) and pre-MS stages. Due to the density and effective 
charge dependence of the opacity enhancement for ModK8, 
during the MS stage the tracks for a Cen A are much closer 
to each other than those of a Cen B. For a Cen A, both 
the standard and improved models are in good agreement 
with the observed results. However, the tracks near the zero- 
age MS (ZAMS) phase and in the greater part of the MS 
phase are slightly different. The effective temperature dif- 
ference between the tracks during the MS (except near the 
terminal-age MS) phase is about 40 K. For a Cen B, how- 
ever, the difference is about 150 K. The improved model of 
a Cen B is in good agreement with the observed position. 
The luminosity difference between the ZAMS of the tracks is 
about 0.04 L and reaches 0.06 Lq for the present models. 
For a given value of luminosity, the radius of the improved 
model of a Cen B is 4 per cent greater than that of the stan- 
dard model. In the following subsection, we consider stellar 
models with mass lower than that of a Cen B. 



5 PRESSURE IONIZATION AND OPACITY 
ENHANCEMENT 

It is very difficult to reveal the physics underlying opacity en- 
hancement. However, it seems that the major role is played 
by the non-ideal effects of pressure (density) ionization. 

Our EOS takes pressure ionization into account by us- 
ing the classical Debye-Hiickel approximation (DHA) (Yddiz 
and Kiziloglu 1997). Detailed analysis of Coulomb inter- 
actions between ion-ion, ion-electron and electron-electron 
shows that the DHA overestimates the Coulomb effect. 
SCVH, for example, consider the ionization of hydrogen at 
low temperatures and increasing densities via a first-order 
'plasma phase transition' and find the existence of a signifi- 
cant amount of neutral hydrogen (5 per cent) at the bottom 
of the convective zone, where the DHA predicts complete 
ionization (see fig. 24 of SCVH). Such a difference is not 
very important for the pressure (fig. 25 of SCVH); the max- 
imum difference between pressures of the two EOS is about 
1 per cent. For opacity, however, a small amount of neutral 
hydrogen plays a significant role. The neutral hydrogen be- 
neath the solar convective zone extends down to the point at 
which log(T) ~ 6.7. This is the region in which the standard 
and the enhanced opacities for the solar models are different 
(see Fig. 5). 

The EOS (Rogers, Swenson & Iglesias 1996) used for 
OPAL opacity tables, however, is different from the MHD 
EOS, and based on an expansion of physical theories (see, 
SCVH and Rogers & Nayfonov 2002). This expansion cannot 
converge in the pressure-ionization regime and consequently 
reduces the sensitivity of OPAL opacity tables. 



4.3 Influence of opacity enhancement in structure 
of 0.8 and 0.9 M models 

In order to test influence of the opacity enhancement on 
the structure of stars cooler than a Cen B, we construct 
0.8 and 0.9 M Q models with ModK8zl9 values of a Cen 
B: X = 0.7126, Z = 0.019 and a = 1.874. In Fig. 8, the 
evolutionary tracks of these models with enhanced (solid 
lines) and standard (dotted line) opacities are plotted in the 
HR diagram. Opacity enhancement seriously decreases the 
ZAMS luminosity; c51ogL = log(i cnh /i st ) is -0.14 for O.8M 
and -0.11 for O.9M , where enh and st represent improved 
and standard models, respectively. Radii of models for a 
given mass is almost the same near the ZAMS. However, 



6 CONCLUSIONS 

Opacity, as one of the main ingredients of stellar modelling, 
directly forms inner and outer layers of the stellar interior. It 
is effective in the nuclear core and stellar envelope because it 
controls energy transport. Therefore, any deviation between 
theoretical and observational results may arise from a pos- 
sible uncertainty in opacity. Considering opacity enhance- 
ment as a function of temperature does not solve the incon- 
sistencies between the models of the Sun and a Cen A and 
B and their seismic and non-seismic constraints. However, 
opacity enhancement treated as a function of density, ion- 
ization degree of heavy elements (oxygen) , and temperature 
(equations 4-5) overcomes these problems. The enhancement 
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below the convective zone, where pressure ionization is ef- 
fective, is about 8k/k = 0.1 for the Sun and 0.17 for a Cen 
B. 

Opacity is influenced by the excitation level and 
ionization degree of ions. These processes are function 
of density and chemical composition as well as temper- 
ature. Therefore, temperature-dependent opacity enhance- 
ment may solve problems for a single star but not for mul- 
tiple stars. The physical conditions for a given temperature 
are different for different stars. The density at log(T)= 6.5, 
for example, is 0.56 g cm~ 3 for a Cen A and 0.65 g cm~ 3 
for a Cen B. The required opacity enhancement at log(T) = 
6.5 for these stars is 7 and 22 per cent, respectively. Increase 
in density at a fixed temperature causes ions to reduce ion- 
ization degree. This process enhances the opacity. 

The evolutionary tracks of cool stars are significantly 
changed by the opacity enhancement (equations 4-5). The 
effective temperature of the enhanced model around ZAMS 
is about 300 K smaller than that of the standard model, be- 
cause luminosity is seriously decreased by the enhancement. 
The effective temperature difference between standard and 
enhanced models with the same luminosity is about 100 K. 
This corresponds to 5 per cent increase in radius. A fur- 
ther test of the opacity enhancement as a function of local 
physical conditions in the radiative regions of cool stars in 
some binary systems is needed. Such an enhancement may 
also be an alternative reason for why variable mixing-length 
parameter is required for the late-type stars of the Hyades. 
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